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Abstract 

We introduce here a simple finite-dimensional feedback control scheme for 
stabilizing solutions of infinite-dimensional dissipative evolution equations, such 
as reaction-diffusion systems, the Navier-Stokes equations and the Kuramoto- 
Sivashinsky equation. The designed feedback control scheme takes advantage 
of the fact that such systems possess finite number of determining parameters 
(degrees of freedom), namely, finite number of determining Fourier modes, deter- 
mining nodes, and determining interpolants and projections. In particular, the 
feedback control scheme uses finitely many of such observables and controllers. 
This observation is of a particular interest since it implies that our approach has 
far more reaching applications, in particular, in data assimilation. Moreover, we 
emphasize that our scheme treats all kinds of the determining projections, as 
well as, the various dissipative equations with one unified approach. However, 
for the sake of simplicity we demonstrate our approach in this paper to a one- 
dimensional reaction-diffusion equation paradigm. 

Keywords. Reaction-diffusion, Navier-Stokes equations, feedback control, data 
assimilation, determining modes, determining nodes, determining volume ele- 
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1 Introduction 

Dissipative dynamical systems, such as the Navier-Stokes equations, the Kuramoto- 
Sivashinsky equation, the complex Ginzburg-Landau equation and various reaction- 
diffusion systems are known to have a finite-dimensional asymptotic (in time) 
behavior (see, e.g., [5], [6], [8], [15], [23], [29], [32], [M], and references therein). 
This is evident due to the fact that such systems possess finite-dimensional global 
attractors ([3], [8], [9], [2U] , [52] , [M] ) , and finite number of determining modes 
(PH> P3> PSIIZZ]), determining nodes ([15], [20], [21], [22], [25], [27], [28]), 
determining volume elements (|22j,|26j) and other finite number of determining 
parameters (degrees of freedom) such as finite elements and other interpolation 
polynomials (|5j,|6j,|21j.) Moreover, some of these systems, which enjoy the prop- 
erty of separation of spatial scales, are also known to have a finite dimensional 
inertial manifolds (see, e.g., [8], [9], [18], pjJ], [31], and references therein). That 
is, in the presence of separation of spatial scales the long-term dynamics of such 
a system is equivalent to that of a finite system of ordinary differential equations. 

There has been some interesting work on reduction methods, with applica- 
tions focused on scientific computing and feedback control theory, taking advan- 
tage of the finite-dimensional asymptotic behavior of these dissipative dynamical 
systems (see, e.g., [1] , [10] , [H] , [21], [33] and references therein) . However, there 
has been very little rigorous analytical work, in particular in the context of 
feedback control theory, justifying these applications. In the case of separation 
of spatial scales, and hence the existence of inertial manifolds, the authors of 
|30j and |31] provide an example of finite-dimensional feedback control (lumped 
feedback control) that drives the dynamics of one-dimensional reaction-diffusion 
system to an a priori specified finite-dimensional dynamics. It is worth stressing 
again that in the case of inertial manifold the dynamics of the underlying evo- 
lution equation is equivalent to that of an ordinary different equations to begin 
with. However, the main challenge is in being able to provide a representation of 
this ODE system in the relevant parameters dictated by the applications. In [22] 
and [6] the authors have shown that if a certain dissipative system has separation 
of scales, and hence an inertial manifold, then such a manifold can be parameter- 
ized by any set of adequate parameters, e.g. Fourier modes, nodal values, local 
volume averages, etc... In the above mentioned work of [30] and [31] the authors 
employed such an equivalence in the parameterization of the inertial manifolds 
to show their results. 

In this paper we propose a new feedback control for controlling general dis- 
sipative evolution equations using any of the determining systems of parameters 
(modes, nodes, volume elements, etc..) without requiring the presence of separa- 
tion in spatial scales, i.e. without assuming the existence of an inertial manifold. 
To fix ideas we demonstrate our idea for a simple reaction diffusion equation, 
the Chafee-Infante equation, which is the real Ginzburg-Landau equation. It 
is worth mentioning, however, that this new idea has a far more reaching ar- 
eas of applications, other than feedback control, such as in data assimilations 
for weather prediction [2J. In addition, one can use this approach to show that 
the long time-time dynamics of the underlying dissipative evolution equation, 
such as the two-dimensional Navier-Stokes equations, can be imbedded in an 
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v Uxx — olu + u 3 = (1) 



infinite-dimensional dynamical system that is induced by an ordinary differential 
equations, named determining form, which is governed by a globally Lipschitz 
vector field, cf. [H], p3] and pi]. 

In this paper we will use the Chafee-Infante reaction-diffusion equation 

du 

u x (0) = u x (L) = (2) 

for a > 0, large enough, as a paradigm to fix ideas and to use the notions 
of finite number of determining modes, nodes and volume elements to design 
feedback control to stabilize the v(x) = unstable steady state solution of ([!])- 
([2]). Indeed, by linearizing equation ([I]) about v = one obtains the linear 
equation 

££ - w xx - av = (3) 
v«(0) = v x {L) = 

Solving the equation ([1]) with initial value vo(x) = A k cos(^7r), with A k € IR, 
and writting a solution of the form v(x,i) = a k {t) cos(^7r) one obtains 

dk + va k C^-) 2 -aa k = 0, (4) 

whose solution is 

a k {t) = A k e^t?-*)\ (5) 

Therefore, for a > 0, large enough, all the low wave numbers k 2 < ^l^l. are un- 
stable. Consequently, the dimension of the unstable manifold of v = behaves 

like J^- (see, for instance, [3], [23] and [Mj for a similar analysis). 
The aim of this paper is to design a feedback control that stabilizes v = 0, for 
example, either by observing the values of the solutions at certain nodal points, 
local averages of the solutions in subintervals of [0,L], or by observing finitely 
many of their Fourier modes. Based on the above discussion, a naive analysis 

would suggest that one would need about J feedback controllers to stabilize 
v = 0. 

In this paper we will give a rigorous justification to this assertion. First, we 
demonstrate our result for the case of local averages, which is the most straight- 
forward approach. Later, we present a more general abstract result, that unifies 
our approach, utilizing all sorts of approximate interpolating polynomials (in- 
terpolants), observables and controllers, and show that this abstract approach 
applies to the Fourier modes, local volume (i. e. local averages) and nodal value 
as particular examples. It is worth mentioning that the same feedback control 
scheme can be used to stabilize any other time-dependent solution of (p])-©. The 
details of the proof are similar to the ones presented here for stabilizing the zero 
solution; thus, for the sake of simplicity they will not be provided. Furthermore, 
similar scheme can be also implemented for feedback control of other nonlin- 
ear dissipative dynamical systems, such as the two-dimensional Navier-Stokes 
equations, the Kuramoto-Sivashinsky equation and reaction-diffusion systems. 
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2 Finite volume Elements Feedback control 

To fix ideas we propose the following feedback control system for {I])-® in order 
to stabilize the steady state solution v = 0, 

n N 

~qI~ v u ™ -otu + u = -fj, } ;u h x Jk (x) (6) 

k=l 

u x (0) = u x (L) = 0, (7) 
where Jk = [(k — k-j^] , for k = 1, . . . , N, Xj ( x ) is the characteristic func- 
tion of the interval Jk, and 

If N f 

¥k = T-j-r / <p(x) dx = — ip(x) dx. 

Wk\ Jj k L Jj k 

Here, the local averages of the solution, Uk, for k = 1, N, are the observables, 
and they are also used as the feedback controllers in ©. It is easy to observe 
v = is also a steady state solution for ([6])- ([7]). 
For ip G H l (\Q,L}) we define 

IMli 1 := J o P 2 (x)dx + J^ (fl(x)dx. (8) 

Before showing that ([6j)-([TJ) globally stabilizes the steady state v = 0, one has 
to prove first the global existence and uniqueness of the feedback system ©-([7]). 
In section [H we will show in Theorem 14.11 a result concerning global existence 
and uniqueness for a general family of finite-dimensional feedback control that 
includes system (|6])-(j7|) as a particular case. Therefore, we will postpone this 
task of proving the global existence and uniqueness until section |33 and we only 
show here the global stability of v = 0. This is in order to fix ideas and to 
demonstrate our general approach. 

Next, assuming the global existence and uniqueness of ©-d?]), we will show 
that every solution u of ©-([T]) tends to zero, as t — > oo, under specific explicit 
assumptions on N,i/,ct,L and /i (see Theorem 12.11 for details). First, we state 
the following proposition which is needed for our result. It is worth mentioning 
that similar propositions were introduced and proved in [7], [22], [25], [26] and 
[27] (see also [30J and [31]). We adapt here similar ideas from [7] for our proof. 

Proposition 2.1. Let ip € i7 1 ([0, L]) then 

N 

MO - Yl ^ kX ,jS'^ L2 - ^H^ll^ 2 - hW^Wm, (9) 



k= 



where h = -j^ . Moreover, 



\v\\h^(J^) (7 2 M + Ileitis), (io) 



where 



N 



7 2 (^) = E 



a- 



k=l 
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N ^ L / N \ 2 

fc=i fc V fe=i fc / 

L / N N \ 2 

= J q [v(x)Yx Jk (x)-Yv k X Jk (x)\ dx, 

\ fc — 1 1 / 

where in the last equality we used the fact that £)fcLi X (a;) = 1, almost every- 
where. Therefore, 

fc == 1 \ fc — 1 / \ / — X / 

= / E (^( a 

■ /o W =i 



fa) - ^z)x 7 fa)x 7 (x)dx. 

Jk Jl 



Since x ( x ) X ( x ) — X ( x )$kli it follows from the above that 

Jl Jk Jk 



M-)-Y,^x k (-)\\h = / Efafa)-^) *jS x)dx 
k=i J ° \k=i j k 

- r 

E / (^)-^fc) 2 dx. (11) 
t~, J Jh 



k=l 

By vertue of Poincare inequality we have 

2 



jf fafa) -^ fc ) 2 < (A) jf (<//(*)) 2 dx. (12) 



'J fc JJ k 
Thus, (HU) and (O imply 

M-)-J2 1 PkX J {-)\\h < f E/ (^'fa)) 2 ^ 



y/(x)) dx, (13) 



which proves inequality Q in the Proposition 12.11 

Next, we prove inequaliy (jlOp . From the Poincare inequality (|12p we have 



jT V\x)dx-T? k h<(^j (v'(x)f dx. (14) 

Thus, by summing over fc = 1, . . . , N, in the above inequality we conclude in- 
equality (fl~0|) of the Proposition 12.11 
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Theorem 2.1. Let N and \i be large enough such that [i>v> (■£-) max{a,fj>} 
where a > and h = ^. Then \\u(t)\\ L 2 tends to zero, as t —> oo, for every 
solution u(t) of 

Proof. Taking the L 2 inner product of equation ([6]) with u, and integrating 
by part, gives 

1 d | - 1 1 o || || 2 1 1 1 1 2 1 1 1 1 4 M V ^ L o U L 2 / \ 

2 ^NIl* + " KIL* - allull^ + H^ = -_^_ = --- 7 („), 

i=i 

and a result we have 

1 d „ „n „9 / /l 



"Ilia- 



Using (I10p . from Proposition I2.1[ and the assumption fj,h > v we have 

"Kllla+M 7 2 (u) = v(\\u x \\ 2 L2 +j 2 (u)) + (v (^j-u)-y 2 (u) > 

(16) 

Substituting (fT6j) in (fT5j) we obtain 

1 d .. .., .^47r 2 ... 2 
2^ IHIi2+( ^- a)IHIi2 -°- 

Therefore, by virtue of Gronwall's inequality and the assumption that v > a 
one obtains 

\\u(t)\\h < e-W^) 2 -«)*||n ( 0)||l 2 ; 
and the Theorem follows. 



Remark 2.1 

It is worth mentioning that if we choose \x = 0(a) then the assumption N > 
in Theorem 12.11 is consistant with the fact that the dimension of the 

unstable manifold about v = is of order of y ^r L - Moreover, in Theorem 15.11 
we give a different and more general proof that illustrate this point further. 



3 Approximate interpolant feedback controllers 

In this section we will consider a general linear map Ih '■ -fT 1 ([0, L\) — > L 2 ([0, L\) 
which is an approximate interpolant of order h of the inclusion map i : H 1 =— )■ L 2 , 
that satisfies the estimate 

11^-4(^)1^2 < ch\\f\\ H i, (17) 

for every tp € i? 1 ([0,L]). The last inequality is a version of the well-known 
Bramble-Hilbert inequality, that usually appears in the context of finite elements 
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[3]. We propose here to consider the following general feedback system of the 
form 

du o 

-q^-vUxx- (xu + u =-nI h (u) (18) 
u x (0) = u x (L) = (19) 

in order to stablize v = 0. Here one can think of Ih(u) as the observables and 
controllers that will be used to stabilize our system. 

Before we state and prove our general theorems concerning system (118p - (ll9p . we 
will give some examples of the approximate interpolant Ih{y) which satisfy the 
approximation property (I17p , In particular, we are interested in approximate 
interpolant 1^ of finite rank, and whose rank is of the order 0(l/h). 

3.1 Examples of finite rank approximate interpolants 

3.1.1 Finite volume elements 

Using the notation of section [2] we consider the approximate interpolant 

N 

4(¥>)=J>kX (z), (20) 
j=i 

that uses local spatial averages (finite volume elements) for approximating the 
local values of the underlying function. We observe that the approximate inter- 
polant Ih(tp) that is introduced in (1201) and implemented in (fTKj) is exactly the 
same one discussed in detail in section (2) In particular, one can easily see that 
approximating inequality (|17h holds in this case, thanks to Proposition 12.11 

3.1.2 Approximate Interpolants based on nodal values 

In this example we consider the approximate interpolant 

N 

h(v) = J2v(xk)x Jk (x), (21) 

k=l 

where Jk and \ are as hi section El and the points Xk € Jk, for k = 1,2, • • • ,N 

Jk 

are arbitrary. Next, we show that approximate interpolant given in (12 1 [) satisfied 
the approximation property ()1T|) . Here again we adopt ideas from [7] to prove 
the next proposition. 

Proposition 3.1. For every ip S // 1 ([0, L]) 

N 

k=l 

Proof. 



~ Yl V^Xj = / fix) - Yl ^( x k)Xj (x) 

k=l Jo \ k =i k J 



2 

dx, 
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N 

and since X h ( x ) = 1 a - e -> it follows that 
k=i k 



N L / N \ 2 

llv(0 - X] V5(^fc)x jft (-)l!i2 = / wZ&ix) ~ <P{xk))x Jk {x))\ dx. 
k=i Jo \k=i k / 

As in the proof of the Proposition l2.lt we observe that \ ( x ) X ( x ) = X ( x ) $ki 

Jk Jl Jk 



and then we obtain 

N . L N 



¥>(•) ~E v(xk)Xj (-)\\h = [ 5^(^0*0 - ^(^)) 2 X Jfe (x)dx 

= E / ~ ^(^/c)) 2 

fc=l ^ J fe 

= E / ( r<p'(y)dy) dx 

k = l ^ ' Jk ^ x k / 

^ E / (7 i^(y)i^) dx 



(22) 



Using Cauchy-Schwarz inequality, we get 

2V AT 



v(.)-J2<p(x k ) Xjk (.)\\h < X> 2 / k(v)| a «fo, 



fe=i fc=i 

= h 2 y x \\ 2 L2 ; (23) 
which concludes the proof of Proposition I3.ll 

In view of (JTJ) , (|18p and (|2ip we propose the following feedback controller for 
stabilizing v = 

^ - vu xx - au + u 3 = -^E n ( x fc) X Jk ( x ) (24) 

k=i k 
u x (0) = u x {L) = 0, (25) 

which is a special case of (fTHj) . 

3.1.3 Projection onto Fourier modes as an approximate inter- 
polant 

Here, we consider the following projection onto the first N Fourier modes as an 
example of an approximate interpolant; 

h(<4>) = E & k COS (— r-)> h = jy » ( 26 ) 
fc=l 
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where the Fourier coefficients are given by 



2 f L ( \ f* kx \A 

LPk = L I cos(—^-)dx. 



Next, we observe that inequality (|17p holds for the approximate interpolant given 
in (BSD. 



Proposition 3.2. Let ip 6 L,L]) be an even function, i.e. (p(—x) = <p(x). 

Then 

HV 9 ^) ~Z2 ^ k COS ( — 7~ J Hi 2 ([0,L]) < c/i ll^||L2([0,L])- (27) 

k=l ^ ' 

Proof. The proof of this proposition is a simple exercise in Fourier series. 
Thus it will be omitted. 



4 Global existence and uniqueness for the 
closed-loop system 

In this section we establish the global existence and uniqueness for the general 
feedback system introduced in (fT8l) - (fT9l) . This will be accomplished under the 
assumption (fT7|) and the condition 

u>iic 2 h 2 . (28) 

To this end one uses the standard Galerkin approximation procedure based on 
the eigenfunctions of the Laplacian, subject to the Neumann boundary condition, 
i.e., cos( 2L ^) for k = 1, 2.... We will omite the details of this standard procedure 
and provide only the formal a-priori estimates (see, e.g., [34 J). These estimates 
can be obtained rigorously through the Galerkin procedure, by passing to the 
limit while using the relevant compactness theorems. 

Let us now establish the aformentioned formal a-priori bounds for the solution 
which are essential for guaranteeing globlal existence and uniqueness. 
System (fl8l) - (fT9l) can be rewritten as 



— - v u xx + j^u - (a + -j-p) u = -u 3 - nh{u) (29) 

u x {0) = u x (L) = 0. (30) 

Taking the L 2 - inner product of (|29]) with u, integrating by parts and using the 
Neumann boundary conditions, we obtain 

1 d f L 2 j f L 2j v[ L 2j f L 4j 

- — / u ax + v u„ ax + —7? u ax = — u ax 



2 dt J i) ,/n f 2 



o Jo ^ Jo Jo 

L r L 



+(ct + Yo) / u 2 dx — n I Ih(u)udx. 
L Jo Jo 
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Writing 

I h {u) u = I h {u) (u - I h (u)) + (I h (u)) 2 
and applying the Cauchy-Schwarz and Young's inequalities, we get 

u 2 dx + v I 
io Jo 



1 d f L o , f L 2 , " f L 2 , 1 f L 4 , / v ,2 L 

— — / u dx + v I u„ dx + — 7 / u dx < — / u dx + (a + — rr - 

2 (ft 7o io ^ 2 Jo ~ 2 y L 2 



- li \I h {u)\ 2 dx + n (J \I h (u)\ 2 dxY \u - I h (u)\ 2 dx" 2 

Using the Young inequality once again we reach 
1 d f L o , /" L 9 . ^ /" L o , . 1 f L A , f , 2 L 



2dt 



I u 2 dx + v I u 2 dx + Y2 I u2 dx<—-f u 4 dx + (a + —^ 
Jo Jo L Jo 2 Jo L 



|4(«)| <ix + - h{u)\\ L 2 dx. 

£ .In ^ 



Using (|17|) . and the definition of the i? 1 -norm given in ([5]) we obtain 

1 d f L 2 , f L 2 , v f L 2 , ^ 1 /' i 4 , , ^ 
/ u dx + v u x dx H =■ / u dx < — / u ax + (a H =■ ) — 

2 Jo Jo L Jo 2 Jo L 2 



u dx + / u^, dx I . 
Jo Jo J 

Thanks to the assumption (j28|) we conclude 

I d f L 2 1 V ( f L 2 , 1 /- £ 2 , \ / " n2^ 

2dt/ n " dX+ 2l/ n ^ + £2 / - ^ <(a + I ,) - 



2 

Thus, one has 



1 f L V f L 

u A dx / \Ih{u)\ 2 dx. 



d " ■ 1 2 , f\\ II 2 1 II ||2 



- \\uf L2 + i/ (H^ll^ + jj- 2 \\u\\ z L2 ) < (a + ) 2 L. (31) 

Therefore, by dropping the ||it x ||^2 term from the left-hand side of (|3ip and 
applying Gronwall's inequality we have 



Ht)\\h < ||u(0)||| 2 + (« + — ) 2 — (1-eT^) =: (32) 



V 2 ^ 3 / 

i 2 ^ t; 

Observe that by integrating equation (j3T|) over the interval [t,t + 1], for t > 0, 
and using ([32]) one has 



t+i i v 

(\\Ms)\\h + j2\\ U ( S )\\h) ds ^ {a + j^) 2 L+\\u(c )U2 

< (a + ^) 2 L + K (t). (33) 
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Next, we find estimate for ||tiir|| 2 . Multiplying equation (|18p by —u xx and inte- 
grating by parts, using the Neumann boundary conditions (fl"9l) . we obtain 



1 d 



f'L f'L fL fL fL 

I u 2 x dx + v I u 2 xx dx = a I u 2 dx + u 3 u xx dx + \i I Ih{u)u xx dx 
Jo Jo Jo Jo Jo 

= a u x dx — 3 / u 2 u x dx + fi u u xx dx 
Jo Jo Jo 

+H \ (Ifc(u) - u) u xx dx. 
Jo 



2dt 1 U - dx + 



Integrating by part, using (|19p. and using the Cauchy-Schwarz and Young in- 
equalities we reach 

-— f u 2 dx + v f u 2 „dx < a f u 2 dx — 3 f u 2 u 2 dx — u [ u 2 dx 
2dtJ x 7 Jo Jo Jo 



ytX 2 .. , . ||2 V 2 



By the assumption (|17p we get 

1 d f L 2 v [ L 2 f L 2 f L 2 2 f L 2 

— — I u x dx + — / u xx dx < a u x dx — 3 / u u x dx — ^ u x dx 

2 dt J Q 2 J Q Jq J Jq 

n 2 

i " 2 i 2 II 1 1 2 

+- C h H H1 . 

Dropping the negative terms in the right-hand side and using the definition of 
if 1 -norm given in ([5]), will lead to 



1 d f L 2 , v f L 2 , , „2 /i 2 c 2 /i 2 / 1 „ ll2 . Il2 

2di / u * dx+ 2 Uxx ~ a ^ Ux ^^ + 2zv 1 l2"H n ^ 2 + ll^^ 2 

^ 2 c 2 /l 2 , 



Therefore, from the above and (1341) we have 



^N|| 2 + |K||£ 2 ). (34) 



|||^|| 2 <2(a + ^|^) (± H 2 2 + IKII 2 ,, ) . (35) 



Since Ko(^)i which is given in ([32]) . is a bounded uniformly, for all t > 0, then by 
virtue of Gronwall's inequality one concludes from (|35p that ||ua;(t)|| 2 2 remains 
bounded on every finite interval [0, T] provided u(Q) € H . However, this estimate 
will grow exponentially in time. Our immediate goal will be to provide a uniform 
bound for ||ua;(t)|| 2 2, for all t > 0. 

From the above discussion, following equation (I35p . we know that one has a 
uniform bound over the time interval [0, 1] . Therefore, we will be focusing on the 
time interval [1, oo). Now, let t + 1 > s > t > 0. We integrate (|35[) over the 
interval [s,t + 1], and use (j33j) to obtain 



2 „2 c 2 ^2 

Mt + l)||| 2 < K(s)|||2 + -(a + ^^) (a + -^)L + K (t) 



Z 2 "- 
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Next, we integrate the above inequality with respect to s over the interval [t, t+1], 
and use (f33l) . to conclude 



K(i + l)||| 2 <- [(a + L + K (t) 



1 + 2(a + 



fi 2 c 2 h\ 
2u ' 



(36) 



for all t > 0. Therefore the if 1 — norm of u is bounded uniformly for all t > 0. 
Next, we show the continuous dependence of the solutions of (|18p on the initial 
data and the uniqueness, provided the assumptions (fT7|) and (f28|) hold. Indeed, 
let u\,U2 be two solutions and w = ui — U2 of (fT8|) . From (JTHJ) we find that 

9io 



at 



fWxa; ~ aw = ul - uf - fl h(w). 



Multiplying by w and integrating with respect to x over [0, L] we get 



1 d 
2~dt 



f 2 f 2 / 2 /" 

/ w dx + v I w x dx = a I w dx — J 
Jo Jo Jo Jo 



2 (Ul + U 2 ) 2 +u\+u\ , 

u; dx 



-fx I Ih(w) w dx 

cL 







< {a. — fi) I w dx + fi \Ih{w) — w\\w\dx. 
Jo Jo 

A straightforward computation, using the Cauchy-Schwarz and Young inequali- 
ties and assumption (fT7|) . yields 



1 d f f f 

-— I w 2 dx + v \ w 2 dx < (a — fi) I w 2 dx + fi\\Ih(w) — w\\ L 2 \\w\ 

2 at J J J 



L 2 



< (a- fx) \\w\\ 2 L 2 dx + ~ \\I h (w) - w\\ 2 L 2 + ^\\w\\ 2 L 2 

M\ ii.. .m2 , M 2 j,2i 



< (a - -) ||w|| L2 + -c A IMIffi 



Using (|8|), the definition of the i/ 1 -norm, we reach 



Id. 

/ w dx + v 

2dt 



w 2 dx 



< ( a - ~ ) 1 1 it? | j 2,2 + ~<r/t 



li lie 2 h 2 



V „2 l2 



L 2 



2 

12 + Ikxiiis 



llz i r" z uz II 



L 2 ' 



By assumption (|28p the above implies 
1 d 



- — IMIl 2 + - < ( a - — + 



2 dV 



x\\L 2 



fl v c 2 h 2 



2L 2 



\w\ 



L 2 ' 



Therefore 



Id. 2 V n2 
2dt M ^ + 2 lML2 - ia 



fi vc 2 h 2 



2 + 2L 2 



w 



|| 2 =: f3\\w\\ 2 L2 , 



and by Gronwall's inequality we have 

IKt)lli 2 <e^|K0)||| 2 . 



(37) 
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Thus, if w(0) = then ||u)(t)||xa = 0. Moreover, inequaliy (|3"7|) implies the 
continuous dependence of the solutions of (fT8j) - (fT9j) on the initial data. In con- 
clusion, from the above, and in particular thanks to f|32j) and (|36p . we have the 
following theorem: 

Theorem 4.1. Let fj,, v and h be positive parameters satisfying assumption &28\) : 
and that 1^ satisfies JT7] ]. Suppose T > and uq G -ff 1 ([0, L]), then system 175])- 
[TV\) has a unique solution u G C([0,T], H 1 ) n L 2 ([0, T], if 2 ) which also depends 
continuously on the initial data. Moreover, 

limsup ||u(i)|||2 < Rq ■= (a + -p^f — , 

and 

limsup ||Wa;(*)|lx,2 < -R? := — 
t-s>oo ^ 



1 + 2(q 



^h 2 
2v 



5 Stabilization using as a feedback con- 
trol 

In the previous section we have established the global existence and uniqueness of 
the general feedback system (fT8|) - (fT9j) . under the assumptions (fTTJ) and ([28]) . In 
addition, we gave in section [3] specific examples for the approxiamate interpolant 
Ih satisfying (|17p . Now, we are ready to state and prove our main result con- 
cerning the general stabilizing feedback system (|18l )-( fT9l) . under the assumptions 
([T7|l . (|25 ]1 and ((58jl below. 

Theorem 5.1. Let I h : F x ([0,L]) -> L 2 ([0,L]) 6e a linear map, which is an 
approximate interpolant of order h of the inclusion map i : H 1 ■-)■ L 2 , i/iai 
satisfies the approximation inequaliy (fj7| ). Moreover, assume that \i is large 
enough such that 

r:=/j-(2a + ^)>0 (38) 
and t/iai /i is small enough such that It28\) is satisfied. 

Then, for every uq € H 1 , the global unique solution of [T8\) - [W\) decays exponen- 
tially to zero as described in |^ j[ ), below. 

Proof. We present here a slightly different proof than the one we gave for 
Theorem 12 .!( taking into account the general form of the approximation property 
(|17|) . However, we emphasize that the main idea is similar. This proof will also 
require simpler explicit choices for the parameters fx > and h. 
First, let us observe that assumption (f28l) guarantees the global existence and 
uniqueness of the solutions to system (fT8|) - (fT9|) . because of Theorem 14.11 
We take the L 2 inner product of (I18j) with u to obtain 

1 d 
2dt 



L rL rL rL 

-— I u 2 dx + v \ u 2 dx — a u 2 dx + / u 4 dx 



J o Jo Jo 

f L 

= —fx / Lh{u)udx. (39) 
Jo 
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From property (|17p . and the definition of the H 1 — norm, given in ([8]), we have 
-2 J I h {u)udx< c *h 2 (\\u x f L2 + ^\\u\\l^J -\\u\\ 2 L2 - \\I h {u)\\ 2 L2 . (40) 



Substituting (@DJ) into pJj) to obtain 
1 d 



ii n2 /" 4 j ii n2 

a II u IIl 2 + / u "- x + tt IfIIl 2 
Jo 2 



Thanks to assumption (l28l) the above implies 

1 j, mo ^ ii no ii t / \n2 ii n4 

ML 2 + o IKHi 2 + oll4(«)|| L 2 + IMIx* 



2 eft" 2 

/,, ,.„2j,2\ 

2 



f 2 " " " "^p- ) II^Hl 2 ^ °- 



By assumption ([28]) we get 
1 d 



2 ^ M u||£ 2 + - Huxllla + §P/»(«)ll2» + \\u\\l<. + - |M|£ 2 < 0. 



Thanks to assumption (J38J) and Gronwall's inequality we have 

\\u(t)\\ 2 L2 <e- rt \H0)\\h- (41) 

Remark 5.1 

Let us observe that in order to satisfy assumption ([38]) one can choose fi = 0(a). 

As a result, asssumption ([28]) will hold if we choose N := ^ = O(y^p-), that 
is the number of feedback controllers is comparable to the dimension of the 
unstable manifold about v = 0. This is consistent with our earlier observation 
in the inroduction and in Remark 2.1. 

6 Stabilizing in the i^-norm 

In the previous section we have shown that the feedback system (|18 p - (|19p sta- 
bilzes the steady state solution v = in the L 2 -norm, i.e., ||tt||x,2 — > 0, as t — > oo, 
provided assumptions ([28]) and ([38]) hold. 

Next, we show that we also have ||w(t)||#i — >■ 0, as t — > oo. To this end it is 
enough to show that H^Hl 2 — > 0, as t — > oo. 
Let us rewrite (fT5D-(fT5H as 



1 1 3 

Ut + J2 U ~ VUxx ~ (° + J2^ U + U = ~f J ' I hW (42) 

u x {0) = u x {L) = 0; (43) 
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We take the 1? inner product of (|42p with —u xx . Notice that from the proof 
of global existence and uniqueness of Theorem 14. 1\ in particular from the uni- 
form boundedness of the H 1 — norm and equation (|34p we can show that u € 
L 2 ([0, T],H 2 ) , and consequently |ei 2 ([0, T],L 2 ) . Therefore, the above men- 
tioned inner product, between ^ and — u xx makes sense rigorously. Integrating 



by parts, and using the Neumann boundary conditions (|19p we obtain: 

^IWli 2 + J2 W U A ~ (« + ^)IM| 2 

u 3 u xx dx + fi I h (u) u xx dx 



= —3 / u 2 u 2 x dx + [i I (Ih(u) — u) u xx dx 
Jo Jo 

+ n I u u xx dx 
Jo 

= —3 / u 2 u 2 x dx + [i I (Ih(u) — u) u xx dx 



o 



fi u 2 dx 



o 



By Cauchy-Schwarz inequality we have 



Id,, ,,9 ,, ,,9 1 

■\\ux\r L 2 

+ H\\lh(u) - U) \\ L 2 \\U XX \\ L 2. 



2^\\ u *\\h + u \\ u xx\\l2 + j2\\ u x\\h - ( a + fp) IKIli 2 ~ A* IKIli 2 



Applying Young's inequality we obtain 

T^IKHi 2 + v \ U *A ^ (« + J2 - tiWuxW 2 + ^ \\u X x\\ 2 L 2 + ^ - n)||| 2 . 

Using property ()1T|) and the definition of the i/ 1 — norm in (|8|), we have 



Id o Z/ 9 1 ||9 1 

.2 ,,2 „2 



+ [ T2 \\ux\\ L 2 + \\u\\ L2 



H,, „, , i/ „ , l0 1 „ „9 , 1 /i 2 " 2 " 2 



2^ir^il! 2 dx + ^ ll n ^ll! 2 + j^ll^lli 2 ^ + -m) IK Ilia 

Notice that ||u x || 2 2 = — J L uu xx dx < \\u\\ L 2 \\u xx \\ L 2. Thus from above we have 
1 d 2 v 2 1 2 1 h 2 fi 2 c 2 2 1.2 

2^ll^^ 2 + 2 H Uxx ^ 2 + pIMI^ - ' a + i2 + ^; — ^1 II m ^Hl 2 IHIl 2 
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By Young's inequality we have 



Id., Il2 , v.. Il2 1 ,, . . 2 (, 1 h 2 u 2 c 2 , 2 1 h 2 u 2 c 2 ' 

^j t \\u x \\^dx + -\\u xx \\ L2 + ^\\u x \\ L2 < n a +- + — —-»)- + — 2 



2 
|| ? 



which implies 

1 d 

2dt u " XUL ^~ ' 4 '™'^ ' L 2 

h? p? <? \2 1 , h? p? (? 



U X \\ L 2 dx + - 11^11^2 + 72 \Wx\\ L 2 < l\\u\\ L 2, 



where 7 = [{a + 7^2 H /^) 2 ^ + ~2UT r ) ■ Consequently, we have 

1 d 11 ■ 1 2 r, , 2 ||„, ||2 / \\„.U\\\2 



X \\ L 2 ux -r 72 U u x\\L 2 - l\\ u \ b )\\L 2 - 



2dt n xnLZ ^ ' L 2 
Since 

lim ||u(t)|| 2 = 0, 



t— too 

2 



then, by Gronwall's inequality it is easy to show that ||«a;||^2 - > 0, as t — > 00 
(see, also special Gronwall's type Lemma in [2S]). Notice that the sign of 7 does 
not matter in the above argument. 



7 Nodal observables and feedback controllers 

In this section we propose a different feedback control based on nodal value 
observables and feedback controllers. Assume that the observables are the values 
of the solutions u(Wk), at the points x k G Jk = [(& — 1) j?t k j?], k = 1, N, and 
that the feedback is at some points xj~ G Jk, %k is n °t necessarily the same as 
Xk. That is the measurements are made at x k , while the feedback controllers are 
at Xk, for k = 1, 2..., N. To avoid technical issues that are dealing with boundary 
conditions, we focus here on the periodic boundary condition case. In this case 
the feedback system will read 

Tj^ - vu xx - au + u 3 = -(jl ^2 hu(x k )5{x - x k ), (44) 

k=i 

u(x,t) = u(x + L,t), (45) 

where h = ^; and 5(x — a) G iT~J.([0, -L]), for a G [0,L], and is extended 
periodically such that 

< £(• - a), <p >= ip{a) (46) 

for every ip G i^ er ([0,L]). 

The feedback control proposed in (gl]) -(|4"5j) is different than that of (fT8|) - (fT9|) . 
since the right-hand side in (|44l) is a distribution that belongs to i^J,([0,L]), 
while the right-hand side in (fl~8|) belongs to L 2 ([0, L]). 

In this section we will show that, under similar assumptions to those in Theorem 
15.11 the proposed feedback system ([31]) stabilizes the steady state v = in the 
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L 2 — norm. One should not expect here a stronger statement, as the one stated in 
section El in which the stabilizing is also valid in the H 1 — norm. This is because 
the solutions of (|44j) -(|45j) are weaker than those of (|18p - (|19p . since the right-hand 
side in (|44j) is less regular than its counterpart in (|18j) . 

Below, we will show the formal steps, which demonstrate simultaneously the 
global existence and stability. These formal steps and estimates can be justified 
rigorously by implementing the Galerkin procedure based on the eigenfunction of 
the Laplacian, subject to periodic boundary conditions, with period L (see, e.g., 
[34J. First, let us prove the following Lemma, which is basically the embedding 
of the Holder space of Cz CH 1 (see also [7]). 



Lemma 7.1. Let Xk,Xk G J k = [(k — l)h,kh],k 
N G Z + . Then for every 92 G -ff 1 ([0, L\) we have 



1, .., N, where h = j? , 



N 



\f(xk) - y{xk)\ 2 < h \\<p x \\ 2 L 2, 



k=i 



and 



\<f\\h < 2 



N 



h ^ \u>(xk)\ 2 + h 2 \\ip x \\ L 2 



(47) 



(45 



Proof. We prove inequality (|4"7]) for ip G C 1 ([0, L]), and by the density of 
C 1 C H 1 the result follows for every tp G H . 



\tp(x k ) - ip{x k )\ 2 





r x k 




< 


/ p'(s) ds 










< 


Jk\ [ W'{s) 


\ 2 ds = h 




JJk 





\(p (s)| ds 



■Ik 



By summing the above inequality over k = 1,..,N we conclude (|47p . 
To prove (|48p we observe that for every x G Jk we have 



|^)| < \<p(x k )\+ / \<p'(s)\ds. 

JJk 



Thus 



\v{x)\ 2 < 2 



\v{x k )\ 2 + 



Jk 



\<p'(s)\ ds 



(49) 



and by integrating with respect to x over Jk, and using the Cauchy-Schwarz 
inequality, we obtain 



Jk 



\ip{x)\ 2 dx <2h 



\<p(x k )\ 2 + h / W{s)\ 2 ds 



Jk 



(50) 



Now we conclude (|48p by summing over k = 1, . . . , N. 



Theorem 7.1. Let /j, > 4 a and h is small enough such that v > 2fih 2 . Then 
for every T > 0, and every uq G L 2 er [0, L] system ^44\ ) has a unique solution 

u G C([0, T];L 2 er [0, L]) n L 2 ([0, T]; H x per [0, L])nl 4 ([0, T]; L^ er [0, L\) , 
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and 

^ G^([0,T];F-i). 

Moreover, 

lim \\u(t)\\ L 2 =0. (51) 

t— >oo 

Proof. We take the H~ l action of (I44p on u € H 1 , and use Lemma of 
Lions-Magenes (cf. Chap. III-p.169, [35J), to obtain: 



1 A r L N 

ii n2 ii n2 n nf 2 / 4 t v — _ 

o~j~IMI2 2 + v II^IIl 2 = a \\ u \\2 ~ u dx — fin / u(xk) u{xk 

_■ City 



L N , N 

« ll^llia - / u 4 dx - nh^2 \u(x k )\ 2 + fj,h ^ (u(x k ) - u(x k )) u(x k ) 
Jo k=i fc=i 



r L N N 

< a\\u\\ 2 L 2- u 4 dx - - h \u(x k )\ 2 + - h ^ \u{x k ) - u(x k )\ 2 , 
J ° k=i k=i 



where in the last step we applied the Young's inequality. Next, we apply ([45 
and (|46p to the right-hand side 



\ d .. . . o Q II mo /* 4 t 7 1 1 ii2 7 2 ii ii2 

«|| L 2 + ^ ||n x .|| L 2 < a ||n|| L 2 — / u ax — — n |[u|| L2 + /i n ||itr|| L 2. 



2 eft 
Hence 

^II u IIl2 +( zy -/ i/l2 )ll^llL2 < (a-^)||«l!i2 n 4 (fx. (52) 

Since v > 2/ih 2 and 4a < « we conclude (|52p . thanks to Gronwall's inequality. 
Moreover, by integrating ([52]) over [0, T], we conclude the regularity of the solu- 
tion as stated in the theorem. Next, we prove the uniqueness of solutions. 
Let «i and U2 be any two solutions. Denote by w = u\ — U2- Then w satisfies 

dw N 

— v w xx — a w + (u — 1 + ui U2 + U2) w = — u h w(x k ) 5(x — x k ) (53) 

k=i 

taking the H~ l action on w G H 1 , and using again Lemma of Lions-Magenes 
(cf. Chap. III-p.169, [35]). we obtain 

Id " L N 

11 ' 1 2 1 11 ||2 || ||2 



/ (u 2 +ui U2+U2) w 2 dx—fi h w(x k ) w(a 
Jo ' k=i 



.-, , u w \\l 2 ~^ v ' \\ w x\\L a ~ a IMIl 2 = / (Ui+ui U2+U2) w dx—puh y ^ w{x k )w{x k ) 

(lb 

Since (u 2 + u\ 112 + U2) w 2 dx > Ul ^" 2 w 2 dx > we obtain 



IQ V I 1 i * 1 ^/ — Jo 2 

TV 

Il7/;ll?i -I- ;y — rv ll?/;ll?o < — /;, h 

2dt 



r— IMI^a + v \w x \^ L i - a \\w\\ 2 L 2 < — iih ^ w(x k )w(x k ) 

k=l 
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From here we follow the same steps as in the proof of the stability to obtain 

1 u ii |,o . ? 2 \ 1 1 ii 2 / U\ || mo 

2di" W " L2 + ( v ~t J 'h ) \\w x \\ L 2 <(«--) \\w\\ L 2. 
Since v > 2fj,h 2 and [i > 4 a we conclude ()52|) . thanks to Gronwall's inequality. 

\\w{t)\\l 2 <e^-^%wmh- (54) 

Notice that (|54|) implies the uniqueness of the solutions and their continuous 
dependence on the initia data. 

Remark 7.1 

Here again we observe that by choosing \i = 0{a) then the condition of the 
theorem imply that N := j t = 0(a/^j) which is comparable to the dimension of 
the unstable manifold about v = 0. 
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